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1 Introduction 

Classically, fundamental solutions are defined for systems of linear partial differential 
equations in M”. Specifically, a fundamental solution to the Stokes system (n > 2) 
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with unknowns v : M” —)■ M”, p : —)■ M and data / : —>■ M”, is a tensor-field 
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We make use of the Einstein summation convention and implicitly sum over all repeated indices. 
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where 6ij and (5]Rn denotes the Kronecker delta and delta distribution, respectively. For 
arbitrary / G a solution {v,p) to (1.1) is then given by the componentwise 

convolution 

Q := * /, (1.3) 

which at the outset is well-defined in the sense of distributions. In the specific case 
of the Stokes fundamental solution Fstokes above, L'^-integrability and pointwise decay 
estimates for {v,p) can be established from (1.3). We refer to the standard literature 
such as [3] and [7] for these well-known results. 

The aim of this paper is to identify a fundamental solution to the time-periodic Stokes 
system 

' dtu — Au Vp = / in M” X M, 

< divu = 0 in M” X M, (1.4) 

_ u{x, t) = u{x, t + T) 

with unknowns u : M"' x M —)• M"' and p : M"' x M —)• M corresponding to time-periodic 

data / : M"' x M —>■ M"' with the same period, that is, f{x, t) = f{x, t -|- T). Here T G M 

denotes the (fixed) time-period. Moreover, x G M"' and t G M denotes the spatial and 
time variable, respectively. The main objective is to establish a framework which enables 
us to define and identify a fundamental solution /Irpstokes to (1.4) with the property that 
a solution (u, p) is given by a convolution 



(1.5) 


Having obtained this goal, we shall then examine to which extent regularity such as 
L'^-integrability and pointwise estimates of the solution can be derived from (1.5). 

Since time-periodic data / : M” x M —>■ M”, {x,t) —>■ f{x,t) are non-decaying in t, a 
framework based on classical convolution in M” x M cannot be applied. Instead, we refor¬ 
mulate (1.4) as a system of partial differential equations on the locally compact abelian 
group G := M” x M/TZ. More specifically, we exploit that T-time-periodic functions 
can naturally be identified with mappings on the torus group T := M/7~Z in the time 
variable t. In the setting of the Schwartz-Bruhat space d^{G) and corresponding space 
of tempered distributions y'{G), we can then define a fundamental solution /Ippstokes to 
(1.4) as a tensor-field 
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in the sense of =y’^(G)-distributions. A solution to the time-periodic Stokes system (1.4) 
is then given by (1.5), provided the convolution is taken over the group G. 

The aim in the following is to identify a tensor-field /rpstokes £ (^Q'j(n+i)xn satisfying 
(1.7). We shall describe Trpstokes as a sum of the steady-state Stokes fundamental solution 
Tstokes and a remainder part satisfying remarkably good integrability and pointwise decay 
estimates. It is well-known that the components of the velocity part T® G ^'(M"-)’^^"' 
and pressure part 7 ® G =5^'(M"')"' of Tstokes are functions 


rUx) := 


7f(x) : = 
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if n = 2 , 


if n > 3, 


respectively; see for example [3, IV. 2]. Here, cOn denotes the surface area of the (n — 1)- 
dimensional unit sphere in MT'. Our main theorem reads: 

Theorem 1.1. Let n > 2. There is a fundamental solution Tj^pstokes £ ^'((7 )(”■+!)fg 
the time-periodic Stokes equations (1.4) on the form (1.6) satisfying (1.7) and 

■yXPS = jS 

with T-*- G satisfying 

Vg G (^1, G L^GT^^, (1.10) 

c 

Vr G [1, 00 ) Ve > 0 3C > 0 V|x| > e : ||r“*“(x, •)||ir(-]r) < pp;, (1-11) 

Vg G (1, 00 ) 3C > 0 V/ G y{Gr : \\r^ * /||n/M..(G) < II/IIl.(g), (1-12) 

where T denotes the torus group T := 'RfTTj, lx G the constant 1, (5t G ^'(T) 

the Dirac delta distribution on T, * the convolution on G, and IT^4.g^(7) Sobolev 
space of order 2 in x and order 1 in t. 


( 1 . 8 ) 

(1.9) 


Remark 1.2. We shall briefly demonstrate how the fundamental solution (1.8)-(1.9) can 
be applied in a more classical setting of the time-periodic Stokes equations to obtain a 
representation formula, integrability properties and decay estimates of a solution. The 
time-periodic Stokes equations are typically studied in a function analytical framework 
based on the function space 


Co“e,(M" X M) := {/ G ^“(M" x M) | f{x,t + T) = f{x,t) A / G x [0,71)}, 
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upon which \\f\\q := ||/||l>j(R"x[o,T]) is a norm. Time-periodic Lebesgue and Sobolev 
spaces are dehned as 


X M) := X M)" 

X M) := Co“ x 


2 , 1,9 


E ii«“iii? + E ii«Tii 

| c .|<2 |^|<1 


1 
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It is easy to see that Lper(M’^ x M) and iyper’^(M^ x M) are isometrically isomorphic to 
L'?(G) and respectively. Regarding as a tensor-field in Lper(M’^ x M), we 

obtain by Theorem 1.1 for any sufficiently smooth vector-held /, say / G C'^gj,(M”'xM)’^, 
a solution (u, p) to (1.4) given by u ■.= ui + U 2 with 


ui := 


rs 


T 

^ j f{-^s)ds 




T 


1 


U 2 := / — / r {x-y,t-s) f{y, s) dsdy 


(1.13) 


and p(x, t) := [ 7 ® *]Rn /(•, t)] (x). Properties of ui and p can be derived directly from the 
Stokes fundamental solution (T®, 7 ®), which, given the simple structure of (T®, 7 ®), is el¬ 
ementary and can be found in standard literature such as [3] and [7]. To fully understand 
the structure of a time-periodic solution, it therefore remains to investigate U 2 - For this 
purpose, (1.10)-(1.12) of Theorem 1.1 is useful. For example, (1.12) yields integrability 
U 2 G lTper’'^(M"' X M), and from (1.11) the pointwise decay estimate \u 2 {x,t)\ < C\x\~'^ 
can be derived for large values of x. 

Remark 1.3. Theorem 1.1 implies that F-^ decays faster than T® as |x| —>■ 00 ; both in 
terms of summability (1-10) and pointwise (1.11). This information provides us with a 
valuable insight into the asymptotic structure as |x| —?■ 00 of a time-periodic solution to 
the Stokes equations. More precisely, from the representation formula u = ui + U 2 with 
ui and U 2 given by (1.13), and the fact that T-*- decays faster than T® as |x| —)• 00 , it 
follows that the leading term in an asymptotic expansion of u coincides with the leading 
term in the expansion of ui. Since ui is a solution to a steady-state Stokes problem, it 
is well-known how to identify its leading term. In conclusion. Theorem 1.1 tells us that 
time-periodic solutions to the Stokes equations essentially have the same well-known 
asymptotic structure as |x| —)• 00 as steady-state solutions—a nontrivial fact, which is 
not clear at the outset. 

The Stokes system is a linearization of the nonlinear Navier-Stokes system. A funda¬ 
mental solution to the time-periodic Stokes equations can therefore be used to develop 
a linear theory for the time-periodic Navier-Stokes problem. The study of the time- 
periodic Navier-Stokes equations was initiated by Serrin [ 6 ], Prodi [5], and Yudovich 
[ 8 ]. Since then, a number of papers have appeared based on the techniques proposed 
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by these authors. The methods all have in common that the time-periodic problem is 
investigated in a setting of the corresponding initial-value problem, and time-periodicity 
of a solution only established a posterior. With an appropriate time-periodic linear 
theory, a more direct approach to the time-periodic Navier-Stokes problem can be de¬ 
veloped, which may reveal more information on the solutions. The asymptotic structure 
mentioned in Remark 1.3 is but one example. 

2 Preliminaries 

Points in M"' x M are denoted by (x, t) with x G M"" and t G M. We refer to x as the 
spatial and to t as the time variable. 

We denote by Br := balls in centered at 0. Moreover, we let := Bji\Br 

and B^:=W^\B^ 

For a sufficiently regular function u : M” x M —>■ M, we put diU := dx^u. The differential 
operators A, V and div act only in the spatial variables. For example, div u := ^j'^j 

denotes the divergence of u with respect to the x variables. 

We let G denote the group G := M” x T, with T denoting the torus group T := ’^jTTL. 
G is equipped with the quotient topology and differentiable structure inherited from 
M"" X M via the quotient mapping vr : M” x M —)■ G, TT{x,t) := (x, [t]). Clearly, G is a 
locally compact abelian group with Haar measure given by the product of the Lebesgue 
measure dx on and the (normalized) Haar measure dt on T. We implicitly identify T 
with the interval [0,T), whence the (normalized) Haar measure on T is determined by 

r 

V/gG(T): J fdt:=^l f{t)dt. 

T 0 

We identify the dual group G with M” x Z and denote points in G by (^, k). 

We denote by =5^(G) the Schwartz-Bruhat space of generalized Schwartz functions; 
see [2]. By d^'{G) we denote the corresponding space of tempered distributions. The 
Fourier transform on G and its inverse takes the form 

: ^{G) ^(G), ^G[u]{^,k) ■■= j j u{x,t) dtdx, 

R" T 

: ^(G) ^ ^(G), ^-\w]ix,t) := E /d^, 

respectively, provided the Lebesgue measure d.^ is normalized appropriately. By duality, 
=^G extends to a homeomorphism ^g ^ y'{G) —>■ S^'{G). Observe that 

We denote by <5®", (5 t, the Dirac delta distribution on T and Z, respectively. 
Observe that 5z is a function with 6z{k) = 1 if A: = 0 and 5i{k) = 0 otherwise. Also 
note that #t[1t] = 

Given a tensor F G ^'(G)""^™, we define the convolution of F with vector field / G 
y{G)^ as the vector held F =t= / G y'{G)^ with [F =t= f]i := Tij =t= fj. 
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The L'^(G)-spaces with norm H-Hg are dehned in the usual way via the Haar measure 
dxdt on G. We further introduce the Sobolev space 

:= ll/lb, 1 ,, ■■={T. + E 

^i«i<2 m<i 

where (G) denotes the space of smooth functions of compact support on G. 

We emphasize at this point that a framework based on G is a natural setting for the 
time-period Stokes equations. It it easy to see that lifting by the restriction vr|]Rnx[o,r) of 
the quotient mapping provides us with an equivalence between the time-periodic Stokes 
problem (1.4) and the system 



I dtu — Au + Vp = f in G, 

[ divu = 0 in G. 

An immediate advantage obtained by writing the time-periodic Stokes problem as system 
of equations on G is the ability to then apply the Fourier transform and re-write 
the problem in terms of Fourier symbols. We shall take advantage of this possibility in 
the proof of the main theorem below. 

We use the symbol G for all constants. In particular, G may represent different 
constants in the scope of a proof. 


3 Proof of main theorem 


Proof of Theorem 1.1. Put 


:= 





where I G denotes the identity matrix. Since 


M:G^C, M{f,k) 


1 - dzjk) 

\f\^ + i^k 


(3.1) 


(3.2) 


is bounded, that is, M G L°^{G), the inverse Fourier transform in (3.1) is well-dehned 
as a distribution in y'. Now define F^ps ^tps gg (1.8) and (1.9). It is 
then easy to verify that 7 "^^®) is a solution to (1.7). 

It remains to show (1.10)-(1.12). For this purpose, we introduce for k £ Z \ {0} the 
function 





i f \l 

4 i 27r|x| 


n — 2 



(3.3) 
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where Ha ^ denotes the Hankel function of the first kind, and ^/z the square root of z 
with positive imaginary part. As one readily verifies, is a fundamental solution to 
the Helmholtz equation 

- A + = <5Mn in M". (3.4) 

Clearly, Fgg^ G Moreover, its Fourier transform is given by the function 

[Fssn] iO = |^|2 

From the estimates in Lemma 3.1 below, we see that 


/ 



q 

2 


dx < oo 


for q G (l, By Holder’s inequality and Parseval’s theorem, we thus deduce 

- Sz{k)) • rJ'sR(x)](t)|'^dxdt 


T R" 


< c 




k 1 
ssrJ 


R'l T 


dt I dx 


< C 


E 

AiGZ\{0} 


pk 

o< 


. q 
2 \ 2 


dx < oo. 


It is well-known that the Riesz transform iHfc(/) := • =^R"[/]] is bounded on 

L'?(M"') for all q G (l,oo). Consequently, we obtain fRj o 91j(^.jr^[(l — 6i{k)') ■ Fs^sr]) G 
L^{G) for q G (l, Recalling (3.5), we compute 


o^h-^io iH,] 1 [(1 - 5j^{k)) • r^g^]) = r,j 

and conclude (1.10). 

In order to show (1.11), we further introduce 


r^:R^\ {0} ^ C, 



-^log|x| 
1 

(n — 2 )a;^ 


\2—n 


(n = 2), 
(n > 2), 


which is the fundamental solution to the Laplace equation AFr = ^R" in As one may 
verify directly from the pointwise definitions of Fgg^ and Fr, the convolution integral 


I Fr(x - y) rg\^{y) dy =: Fr * Fg^^ix) (3.6) 
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exists for all x G M"" \ {0}. In fact, the function given by 11 * -Tg^gR belongs to 

and defines a tempered distribution in One may further verify that also the 

second order derivatives of 11 * -^I^sr given by convolution integrals 

didjlr, * Tgy (x) = J d,r,{x - v) d,rl,^{y) dy, (3.7) 


from which it follows that their Fourier transform are functions 

We infer from the expression above that 

ri- = 1 [(1 - 5^{k)) • [5,,dhdh - d,dj\ [r, * Fg^g^]]. 

Employing Hausdorff-Young’s inequality in combination with the pointwise estimate 
from Lemma 3.2 below, we obtain for r G [2, oo) 

OllL-m < ( Z] I ■ i^ij^hdh - didj] [Fl * Fg^(x) \ 

^ fcez\{o} ^ 


which concludes (1.11). 

The convolution F"*- * / can be expressed in terms of a Fourier multiplier 




-1 

G 


M{C,k){l- 




^Gif] 


with M given by (3.2). As already mentioned, M G L°°{G). As one may verify, also 
second order spatial derivatives didjM G L°°{G) and the time derivative dtM G L°°{G) 
are bounded. Based on this information, (1.12) can be established. For the details of 
the argument, we refer the reader to [4, Proof of Theorem 4.8]. □ 

Lemma 3.1. The function defined in (3.3) satisfies 

1 

(3.8) 

properties of Hankel functions sum- 
case n > 2. Employing (3.16) with 


<C|A:|-3|xr5 (3.9) 


E \rLi^)\y<G\x 

k£l.\{0} ^ 

Proof. The estimates are based on the asymptotic 
marized in Lemma 3.3 below. We start with the 
e = 1, we deduce 


VA:GZV|x| > 





H 


( 1 ) 

-1 


.27r 

—i — k ■ \x\ 
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Employing (3.17) with i? = 1, we obtain: 


VA:GZV|x| < \ : 

V 27r 
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( 1 ) 

-1 


.2vr 

—i—k ■ \x\ 


n—2 n—2 

< C \k\ |x| 


It follows that 

feez\{o} 


n — 2 r> ^ 2 — n o ^ 

1 ^ \xf-^\k\^ \x\^-^ 




+ E 1^1'"" 1^1"" 1^1"' ) 


l^l>^l*l 
< C( 1 x 1 “^ • 




(ki) 


|fc|>i ^ 


For |g| < 1 we observe that 




OO (j+1)^ —1 ^ 

E E 


k>l 


j=l k=p 

OO 0 + 1)^ —1 

sE E + 

k=p 

OO 

= Eia+ir“'9' 

j=i 


< 9 E-^ (j + 2)... (j + n - 3) ^ 

i=i 


= 9 


in —2 


E«' 


■ i=i 

= (n-2)!. g(l-g)i-", 


from which we deduce 


fcez\{o} 


|2(1—n) 

a^r • Xr 


(N) 


+ |x|^ e 


< C\x 


. g-vrl^l 


(3.10) 


(3.11) 
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and consequently (3.8) in the case n > 2. In the case n = 2, we employ (3.18) to deduce 

(3.12) 


VA; e Z V|x| < \ : 

V 27r 


H, 


( 1 ) 




< C 


log ( ^/y|A:|2|x| 


It follows in the case n = 2 that^ 

|2 


E I^'sh(^)I E 




log ( ^/y|A:Hx| 




< c 


, I /2vr 1. . 

log ( \j —t2|a;| 




(bi) 


(3.13) 


<C\x\ / |log(s)| sds'X 


■^[o.y£ 


(kl) 


<CN- -Xf.^dxl). 

Estimate (3.9) is still valid in the case n = 2. We can thus proceed as in (3.11) and 
obtain (3.8) also in the case n = 2. □ 


Lemma 3.2. The convolution defined in (3.6) satisfies 


Ve > 0 30 0 V|x| > e : |(9i5j[A * r|gg](x)| < C Ixr*". (3.14) 

\x\ 

Proof. Fix e > 0 and consider some x G MT with |x| > e. Put i? := ^. Let x £ C'q 


T(K;] 


be a “cut-off” function with 


Xir) = < 


0 when 0 < |r| < 

- I 1 - 2 ’ 

1 when 1 < |r| < 3, 
0 when 4 < Irl. 


Define XR : M"' —)• M by Xhiv) •= x{R ^|y|)- We use XR to decompose the integral in 
(3.7) as 

didj [El * -Ts^sr] (a;) = f d^T^ix - v) djPssKiy) XrIv) dy 


B. 


‘AR,R/2 


+ / diri^{x-y)djr^sRiy){'^-XR{y))dy 


Bt, 


+ / diR{x - y) djr^s^{y) (l - Xi?(y)) dy 

B3R 

= : Ii{x) Phix) Phix). 


would like to thank Prof. Toshiaki Hishida for suggesting this estimate to me and thereby improving 
my original proof. 
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Recalling the definition (3.3) of as well as the property (3.15) and the estimate 
(3.16) of the Hankel function, we can estimate for \y\ > R/2: 


\dA{y)\<C\k\^ 


di 


2 — n 

\y\^ 


H 


( 1 ) 

n —2 I \/ “ >~r^ 

2 \ V / 
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2 — n 

\y\^ dj 


HFU-i'^k.\y\ 


I ffcO IL ± It ± It ± » 

< C( |A:|4“4 \y\~2-2 + |fc|4-4 |?/|“2 + 2 ) g 

<C\k\-^\y\-^^+^\ 




Consequently, we obtain: 


|/i(x)|<C' / \x-y\^-^\k\-^\y\-^^+^^dy<C\k\-^R 


1-1 D-r 


B. 


4,R,R/2 


To estimate I 2 , we integrate partially and employ polar coordinates to deduce 
\l 2 {x)\<C j |aj5irL(x-y)| |rg^sj,(y)| + |airL(a:-y)| |rs^sR(y)|^"^d?/ 


Br 


<C J R-^\r,%M\dy 

Br 

/ n—2 2—n 

R-^\k\^ \y\^ 

Br 

R 


H 


( 1 ) 


n —z \ \ ^ /-r 

2 \ V / 


■ \y\ 


dy 


< C / R-”|A:|~r 2 


n-2 „ ( 1 ) 

n-2 


-d^k.r 


dr 


< C I R-^\k\-^s 


rr{l) 
ri n-2 


0 


—I 


2tx k 

r'yw\ 


ds. 


Employing in the case n > 2 estimate (3.17) in combination with (3.16), we obtain 

1 00 


|/2(x)| < CR-^\k\-^ 


y—n I 7^1—1 


J s 2 s ds + J s‘2 s 2 e V dsj < C R 
0 1 

When n = 2, we use estimate (3.18) in combination with (3.16) and obtain also in this 
case 



< CR-^\k\-\ 
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In order to estimate /a, we again integrate partially and utilize (3.16): 


|/3(a;)|<C' / \djdirj^{x-y)\\r^g^{y)\ + \dir^{x-y)\\r^s^{y)\R ^dy 

B3R 

<C j R-^\rl,^{y)\dy 

B3R 


n—2 2—n 




dy 


<C / R-^\k\—\y\- 

B3R 

<C f R-^lkl"^ \y\^ dy 

B3R 

<C j R-'^\k\^ \y\^ [\k\^\y\)~'^ dy <C R-^\k\-^^ <C R- 

B3R 


1-1 


Since |x| = 2R, we conclude (3.14) by collecting the estimates for R, R and R. □ 


Lemma 3.3. Hankel functions are analytic in C\ {0} with 


Vz/GC V2GC\{0}: 

) = Hl^},iz)--^Hi^\z). 

(3.15) 

The Hankel functions satisfy the following estimates: 


Vi^ G C Ve > 0 3(7 > 0 V z > e : 

\Hi^Hz)\ <C\z\-'^ e-^“^ 

(3.16) 

Vi/ G M+ Vi? > 0 3(7 > 0 V|z| < R : 

\Hl^\z)\<C\zr, 

(3.17) 

Vi? > 0 3(7 > 0 V|z| < i? : 

\Hl,^\z)\<C\log{\z\)\. 

(3.18) 


Proof. The recurrence relation (3.15) is a well-know property of various Bessel functions; 
see for example [1, 9.1.27]. We refer to [1, 9.2.3] for the asymptotic behaviour (3.16) 
of hI^\z) as z —)■ oo. See [1, 9.1.9 and 9.1.8] for the asymptotic behaviour (3.17) and 
(3.18) of Hi^\z) as z —)■ 0. □ 
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